Instanton liquid at finite quark density 



S. V. Molodtsov, G. M. Zinovjevt 

State Research Center, Institute of Theoretical and Experimental Physics, 117259, Moscow, RUSSIA 

^ Bogolyubov Institute for Theoretical Physics, 
National Academy of Sciences of Ukraine, UA-03143, Kiev, UKRAINE 



Abstract 

The interaction of light quarks and instanton Hquid in the phase of nonzero values of chiral 
condensate is studied at finite quark/baryon density. Calculating the generating functional in 
■ the tadpole approximation we investigate the behaviour of dynamical quark mass and chiral 

! condensate together with instanton liquid density (which shows insignificant increase) as a 

function of quark/baryon chemical potential. We argue on the noticeable increase of quark 
density at the point of developing the nonzero magnitude of diquark condensate due only to the 
quark interaction with the instanton liquid. 
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^ ■ Since the discovery of instanton solutions in the Yang-Mills theory the functional integral language 
T— I , has become very popular and the most operational instrument of quantum chromodynamics (QCD) 
^ ' phenomenology. This formulation tolerates the nonperturbative identification of the field configura- 
tions dominating the action and manifests the fundamental physics features of QCD distinguishing 
o I it from other complex systems with many degrees of freedom. In an overwhelming number of these 
ft) I theories the quanta exchanged between interacting fermions according the initial Lagrangian lead to 
D > the formulation of the concept of interaction potential and the definition of its concrete form. The 
[ dominating role of topological excitations of the gluon field and quantum oscillations about them in 
[ the QCD physics allowed us to understand and successfully explain new physical effects related, first 
^ ' of all, to the process of chiral symmetry breaking. Moreover, it compels us to have some doubts on 
■ the potential concept (perturbative description) pertinent. 

In the intervening years considerable theoretical efforts were devoted to extend the analytical 
description beyond dilute instanton gas approach |jl|. It was shown that an effective instanton- 
anti-instanton interaction could stabilize the medium of these pseudoparticles (PP) @ at the level 
of pretty reasonable values of major parameters (average sizes and distances between PPs) form- 
ing rather instanton liquid The approach developed has succeeded in making interesting (even 
quantitative) predictions although the lattice QCD studies only managed to put the problem of in- 
vestigating the instanton role on the serious theoretical ground IQ. While much has been learned 
about the QCD phase structure at finite temperature, this structure at nonzero quark/baryon den- 
sities has been significantly less explored because of very serious technical difficulties, mainly, with 
putting the fermions (and respective chemical potential) on the lattice. Nowadays, the possibility of 
direct experimental measurements of instanton-induced effects is even discussed in view of the future 
program for deep-inelastic scattering (small size instantons) at HERA P|) and for ultrarelativistic 
heavy ion collisions at RHIC and LHC [Q. 

This paper is addressed to investigate the light quark interactions with the instanton liquid (IL) 
in the phase of nonzero chiral condensate values at finite quark/baryon densities and continue the 
analysis we began in our preceding paper where a new method was developed to estimate the IL 
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back-reaction on the quark presence. Before it was always considered to be negligible. It is natural 
to start with the generating functional of the IL theory having the factorized form 

Z = Zg ■ Z^ , 

where the factor Zg originates information on the gluon condensate while the fermion component Z^ 
provides the opportunity to study the quark behaviour in the instanton medium and, in particular, 
the chiral and diquark condensates together with their excitations The factor Zg is treated in 
quasi-classical approximation supposing the superposition of pseudo-particle (PP) fields Ajj{x]'~^) 
which are the Euclidean solutions of the Yang-Mills equations called the (anti-)instantons (//) to be 
dominant saturating configurations. Such a solution is characterized by the parameters 'j = p, z,U 
where p denotes the PP size centered at the coordinate z with the colour orientation defined by 
the matrix U in the space of colour group SU{Nc). Let us suppose, for clarity, pseudo-particles 
available in four- volume V and a liquid is topologically neutral Nj = Nj = N/2. Formulating the 
variational maximum principle the following practical estimate has been obtained 

— ^ , 

with the average IL action {S) defined by the additive functional 



{S) = J dz J dp n{p) s{p) , (1) 
and with averaging the action per one instanton 

sip) = Pip) + 51n(Ap) - In^S^TV^ + /^^V/ dp, n(pi)p? , (2) 
weighted with instanton size distribution function 

nip) = C e-'^P^ = C p-5^27Vcg-/3(p)-^^pV7 ^ (3) 

where p"^ = J dp p^ nip)/n = [jj^^ , n = J dp nip) = ^, z/ = b = ^ Nf 

is the number of flavours. The constant C is defined by the variational principle in the selfconsistent 

way and Pip) = = —IhCn^ — 61n(Ap) (A = Aj^g = 0.92Ap.y) with another constant C^^ 

depending on the renormalization scheme, in particular, here C^^ ~ ^2^^ exp ^ ^L68A^) ^ ^ r^^^ 

parameters P = Pip) and P = P + IuCat^ are fixed at the characteristic scale of the average PP size 

p . The constant = ~^ j^'i ^ -^ t^"^ gives information on interaction in the stochastic ensemble of 

PPs. Then the IL parameters extracted, for example, average PP size p and the IL density n agree 
reasonably with their values obtained from the phenomenological studies of the QCD vacuum. 

In the quark determinant Z^ the quark fields are considered being influenced by stochastic PP 
ensemble but the quark backreaction upon the instanton medium is usually ignored, i.e. 



Small values of the packing fraction parameter (np^) allow us to neglect the correlations between PPs 
and one usually deals with the limit ^ oo. The dominant contribution to the action of fermion 
fields mainly comes from the zero modes $(x — z; p) which are the solutions of the Dirac equation 

[iDiAjj) - ipii] = , 
with the chemical potential p in the (anti-)instanton field Ajj, and their Weyl components 
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Here ^/j ^ = P± $, P± = r = |a3|, cr^ = (±zcr, 1), cr is the vector of the Pauh matrices and e 

is the antisymmetric tensor. In particular, the quark determinant at Nf = 1 might be written as 



VR. 



VR. 



exp 



dx (i9 — iyU74)'?/'(x) 



(4) 



where the factor R provides the dimensionless result. The preexponential factors are responsible just 
for the instanton-induced interaction of quarks ^ 
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i j dz dU dpn{p)/n I dxdy ipl nix) {id — ifi)^^jj{x — z) ^jj{y — z) [id — ip)^ ipL,R{y) , 



and the colour orientation averaging is performed by integrating over dU . The designation $ is 
introduced for the conjugated zero mode /i) = — /x). /i^ and similar designations are used 

for the four-vectors spanned on the matrices crj^, i.e. /i^ = /i^crj^, = (0,/i). 

Integrating over the auxiliary parameter A brings the quark determinant to the exponential 
form convenient for the saddle point calculation 



r d\ f 

X / Dip^'Di/j exp 
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(^2^^|4^^(^) y-k-ip + j;^lo{k]fi)jij{k)j . 

The complex function 'jQ{k; fi) of this part of the generating functional is defined by the Fourier 
components of the zero modes 

7o(/c; fi) = {k + ifi)4k + hf^ik; fi)hp{k; //) , $(/c; /x)^" = h^{k; /x)(ctJ)} e^'' , 



hiiki^ k; /i) 



~Ak 



{{k — jJ, — ik4)[{2ki + in)fi + i{k — — ik4)f2 ] + 



+ (A; + /i + iki)[{2ki + - + // + iki)f:^]} 



hi{k4,k;fi) 



4P 



{(2A; - - p - iki)f{^ + {2k + p){k + fx + iki)f^+ 



+ 
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2{k - jj,){k - jJ,- iki) - -{jj, + ik4)[kl + {p - yu)^] 
2{k + fi){k + /i + ik^) + hp + iki)[kl + + pf] 
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where k = \k\ for the spatial components of four- vector, 

h{z^)Ko{z^) - Io{z^)K^{z^] 



zi ' 2 



2± 



and li, Ki {i = 0, 1) are the modified Bessel functions. 

In order not to overload the following formulae with the unnecessary coefficients we introduce the 
dimensionless variables 

^-^4, H-A, m^-p'^'m- 

Then it is easy to get the saddle point equation taking the form 

!!2!_2/V =0, (5) 

here the prime is attributed to differentiating in A. 

It has been shown in Refs. 0, that the backreaction of quarks upon IL may be reproduced 



perturbatively with the small variations of the IL parameters Sn and Sp around their equilibrium 
values n and p. These variations are actually accommodated by the approach if the deformable 
(crumpled) (anti-)instantons of size p being the function of x and z, i.e. p — > p{x, z) are utilized as 
the saturating configurations. If the wave length of excitations is much larger than the characteristic 
(anti-)instanton size p (for example, for vr-meson) the mean action per one instanton gains the extra 
term (looking like kinetic energy) generated by the deformable (anti-)instantons 



{S)^ dz dpn{p) - \'^\ +s{p)\ , (6) 




where n is the kinetic coefficient being derived within the quasi-classical approach. If we strive for 
holding the precision declared in this paper the kinetic coefficient should be fixed on a characteristic 
scale, for example, as k ~ '^(p). Our estimates give the value of a few single instanton actions k ~ c /5 
(with the factor c ~ 1.5 — 6 depending on the ansatz for the saturating configurations). Collecting 

the terms of the second order in deviation from the point of action minimum only '^^^^ = 0, and 

taking approximately 

n\ ~ n\ A ^ 

dp 



s{p) ^ s{p) + -4^ ^^ (7) 



where s^'^^{p) ^ ^ ^^2^ 



= , and the scalar field ip = Sp = p — pc — p — p is the field of deviations 

Pc P^ 

from the equilibrium value oi pc = p (1^2^) — P; it becomes clear the deformation field is 
described by the following Lagrangian density [|): 

2 

2 1 \dz 



f^l +MV 



with the mass gap of phononlike excitations = j-^^ = -^=. M ^ 1.21 A for IL in the 

ft. ^p2 

"quenched" approximation with the parameters fixed as A^^c = 3, c = 4, p A ^ 0.37, (3 ~ 
17.5, n A-^ ^ 0.44 |T§. 

Including the variations of zero modes in the quark determinant results in changing those for 
$j/(x - z,p) ~ $/j(x - z,pc) + ^f}{x - z,pc)Sp{x,z), where ^f}{u,pc) = d^ri{u, p)/dp\p=p^ and 

^) It seems to us the physical meaning of the deformation field is analogous to the phonons of solid state physics. 
It is why we called them phononlike excitations of IL . 
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because of the adiabaticity constraint it is valid 6p{x,z) ~ 6p{z,z) = (f{z). The extra contributions 
of scalar field generate the corrections in the kernels of factors which are treated in the linear 
approximation in (p and it is supposed Pc — P everywhere. As result we have the effective Lagrangian 
with the Yukawa interaction of quarks and colourless scalar field 

JdXZgJ Dij^Di; Dip explN (in'^-U-J^^ ip{-k) 4 [P + M^] ip{k)\ x 

(8) 

X exp < y — — ip^k) 2 TT'^6{k — l){—k — ip + iX'yo{k; p)) + iX'yi{k,l; p) ip{k — I) ip{l)> . 

Let us mention the dimensionless variables as the scalar field ip{k) {nK)^^^'^p^ '■p{k) and its mass 

— i> M were introduced here. Zg = e^"^*^***^^)-'^ denotes the part of the gluon component of the 
generating functional which survives expanding the action per one instanton in small deviation from 
the instanton equilibrium size. The general form of the kernel 71 (/c, /; p) is not of crucial importance 
for us here. We may simplify calculation essentially because of the presence of quark condensate and 
collect the leading contributions to the functional Eq. which are just the tadpole graphs. In this 
particular case the following purturbative scheme (in the deviation from the condensate magnitude) 
is applicable 

t/^V ¥^ = (V'V) ¥^ + {V'V - (V'V)} </5 , (9) 



and for the vertex we then have 
2 

{np'^K, 



7i(A;; p) = - — + (^4 + ip)h^][k6h + {kii + ip)5h^ + [(A;4 + ip)h — kh4\[{k4^ + ip)5h — k5h^}. 



Here we introduced the component h spanned on the unit vector — ^ ^ ^.nd Sh^ are 

produced by h and while the substitutions — > gf where 

gf = ft + 2{ltKt - ifKf) , gf = -ff - ff , 

done. 

In the present paper two corrections are taken into account. One of them is related to the changes 
of dynamical quark mass and chiral condensate due to the interaction with the scalar field but another 
one comes from changing the IL density invoked by the shift of the PP equilibrium size. The former 
is given by the following contribution 



^i'>^f I ^16 li{kJ;p)^^{k'J';p) ^\k)m ^\k'm') {V{k~l)^{k'-l')) 



,dk rdl 
4 aV - li{k;p) ^\k)m / - 7i(^;^) Tr S{1) 0(0) , 



where the conventional definitions are introduced for the convolution of scalar field 

{^{k) ifil)) = 7f'6{k + I) D{k) , Z}(A;) = ^-^^^1^ , 

and for the quark Green function S{k) 

{^\k)ij{l)) = -7iH{k - I) Tr S{k) . 

All the factors surrounding ip'^{k)ip{k) in Eq. (|10D give the extra contribution to the quark dynamical 
mass A7o(A;; p) — > Ar(A;; p) = Xjoi^k; p) + A(57o(A;; p), where 

Sjoik; p) = 7i(A;; p) (-2^) / - 7i(/; f^) Tr S{1) D{0) . (11) 
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Let us introduce now the auxiliary function c(A) in order to emphasize the obvious parameter de- 
pendence of 57, i.e. 

If the guark Green function is known c(A) could be defined by the complete integral equation 

c(A) ^ («A)"= I'^M^;,) ■ (") 

what helps also to calculate S{k). 

It was shown in Ref. the variation of equilibrium PP size generates the IL density variation 

and can be given by the tadpole contribution ip{0) A with A = c(A) which produces an extra 

term in the mean action per one instanton 

{S) = Jdzn!^{s)-{A ^)| . 

The maximum variational principle helps to find the PP average size pA = exp | — 2^^^]^ | and to 
obtain the following equation to calculate the IL density 

Eventually, we come to the equation for the saddle point of generating functional Eq. (|^) in the form 



A 7 vr^ {k + ifiY + X^T\k; fi) ^ ' A 

where the IL density variations are incorporated. The prime here signifies the derivative in A. But 
we need to know another derivative c' to have the complete equation. With the definition of c(A) as 
given by Eq. (|T^) , we obtain 



(1 - \^A{\)) c'(A) = 2A A{\) c(A) + B{\) 
where the functions A{X) and -B(A) are the following integrals 



5(A) = -2A {np^'nY'^ J 
and the factor a (A) reads 



u J [{k + ifi)^ + \^T\k; 

dk 7i(A;) V^{k) 



a{X) = 1 - A 



[{k + l^lf + \^v\k■^i)f ' 

^N, V{u+l/2) c(A) 



The corresponding results for multifiavour approach Nf > 1 could be obtained with the simple 
substitutions We have to introduce the factor Nf into the extra contribution to the dynamical 

mass in Eq. (|To|) and to the factor A in the tadpole contribution. Besides, the logarithm in Eq. 



no^ ( fno^V^^^ 1 \ 

(|Tip should be modified according to the prescription In — ^ In ( ( — ^ j 1 . For clarity we 
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Figure 1: The dynamical quark masses (solid lines) and chiral condensates (dashed lines) as the 
functions of chemical potential /i at A = 280 MeV at Nc = 3, Nf = 2. Both lower curves correspond 
to the calculation with the modified generating functional. 



take the factor R to be equal one (see, for example, Eq. (|^)) everywhere. However, this parameter 
could be treated as the free one. Then the dimensional parameters of this approach are given by the 
corresponding powers of Aqcd, for instance, pA and could be selfconsistently obtained if the saddle 
point was already known. 

In Fig. 1 we expose the calculation results with the positive root of Eq. (p^) 



-4 

np 



+ 



+ 



A r(z/+l/2) 
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for the dynamical quark mass Ar(0,/i) (solid lines) and for the chiral condensate —i{ip'^ilj) (dashed 
lines) as the functions of chemical potential yU at A = 280 MeV and for A'^c = 3 , Nf = 2. Both 
lower curves correspond to the calculation with modified generating functional when the tadpole 
contributions are taken into account. Apparently, the results obtained are well within the scope of 
QCD vacuum phenomenology and, in principle, the perfect fit is achievable with slight variation of 
Aqcd- 

Fig. 2 demonstrates rather insignificant change of the IL density when the quarks are in the phase 
where the chiral condensate develops non-zero values though the quark influence on IL is transformed 
into the small strengthening of gluon condensate manifesting itself in the supplementary attraction 
available in the system of quarks and gluons. These features fully agree with the corresponding 
results of Ref. where the IL behaviour was studied in the so-called cocktail model (with (anti- 
)instanton molecula admixed). In that model the IL density in the phase of broken chiral symmetry 
is almost constant up to the critical values of chemical potential. 

The results of studying the saddle point parameter A (which is proportional to the free energy 
within the precision of one loop approximation) at A"c = 3, Nf = 2 are shown in Fig. 3 where the 
upper curve was calculated without corrections coming from the quark feedback taken into account. 
In fact, these results allow us to make one substantial qualitative prediction for the behaviour of 
critical chemical potential pc while transiting to the colour superconductivity phase. Indeed, keeping 
in mind the results of Refs. |]T3[ we should foresee the point pc shifting perceptibly to the larger 
values because the parameter A for the diquark phase at precritical values of p is always disposed 
above than for the corresponding curves of broken chiral symmetry phase. It means the crossing 
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Figure 2: The IL density in the phase of nonzero chiral condensate values a,t Nc — 3, Nf — 2 as 
function of chemical potential /i. 

point with the lower curve is always situated further along the ^-axis than the crossing point with 
the upper curve. Such a situation signals the density of quark matter in the critical interval could 
be comparable (or even larger) than normal nuclear density. In the meantime, we should remember 
that the IL model itself should obey some hmitation for the chemical potential values. It is just 
the region where the transition to perturbative quark-gluon matter phase takes place. With quark 
matter density increasing the average interquark distance may be so small that the interquark gluon 
('Coulomb') field strengths become comparable or even exceed the instanton ones. Then it is invalid 
to consider the (anti-)instanton ensemble as the saturating configuration and the IL approach is 
made irrelevant in this region. 
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Figure 3: Saddle point A as function of n in the phase of broken chiral symmetry at Nc — S, Nf — 2. 
The upper curve corresponds to the solution with no perturbation of IL. 
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